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Abstract 

Endemic species on islands are considered at risk of extinction for several reasons, 
including limited dispersal abilities, small population sizes, and low genetic diver- 
sity. We used mitochondrial DNA (D-Loop) and 17 microsatellite loci to explore 
the evolutionary relationship between an endemic anemonefish, Amphiprion mccul- 
lochi (restricted to isolated locations in subtropical eastern Australia) and its more 
widespread sister species, A. akindynos. A mitochondrial DNA (mtDNA) phylogram 
showed reciprocal monophyly was lacking for the two species, with two supported 
groups, each containing representatives of both species, but no shared haplotypes 
and up to 12 species, but not location-specific management units (MUs) . Population 
genetic analyses suggested evolutionary connectivity among samples of each species 
(mtDNA), while ecological connectivity was only evident among populations of the 
endemic, A. mccullochi. This suggests higher dispersal between endemic anemone- 
fish populations at both evolutionary and ecological timeframes, despite separation 
by hundreds of kilometers. The complex mtDNA structure results from historical 
hybridization and introgression in the evolutionary past of these species, vali- 
dated by msat analyses (NEWHYBRIDS, STRUCTURE, and DAPC). Both species 
had high genetic diversities (mtDNA h > 0.90, n = 4.0%; msat genetic diversity, 
gd > 0.670). While high gd and connectivity reduce extinction risk, identifying and 
protecting populations implicated in generating reticulate structure among these 
species should be a conservation priority. 



Introduction 

Remote islands often contain a disproportionate number of 
endemic species (Ceballos and Brown 1995; Randall 1998; 
Gillespie et al. 2008) and genetically distinct populations of 
species with broad geographic ranges (Avise 1992; Slatkin 
1993). The uniqueness of island communities makes them 
a high conservation priority (Gillespie and Roderick 2002), 
even more so given the extremely high rates of local and global 
extinctions of species inhabiting islands (Steadman 1995). 
Island endemics face an elevated risk of extinction because 
they often have vulnerable traits, such as, limited dispersal 
abilities, small population sizes, and low genetic diversity 
(gd) due to genetic drift and inbreeding (Frankham 1997, 
1998). A sound knowledge of the evolutionary history and 



adaptive capacity of endemic species or isolated populations 
is required in order to understand what increases their risk of 
extinction, thereby enabling the development of appropriate 
conservation measures (Faith 1992; Moritz 2002). 

Throughout the world's tropical oceans, coral reef organ- 
isms are distributed on islands and reefs that represent vary- 
ing degrees of isolation. For coral reef fish, it is the most iso- 
lated locations that contain the highest levels of endemism 
(Lessios et al. 2001; Jones et al. 2002; Allen 2008). Isolated 
islands also support genetically differentiated populations of 
some widespread reef fish species (Muss et al. 2001; Planes 
and Fauvelot 2002; Winters et al. 2010). Islands communi- 
ties are also characterized by a high proportion of vagrants, 
which can result in a high incidence of hybridization due to a 
scarcity of conspecific mates (Hobbs et al. 2009b ). Reef fish on 
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isolated islands also appear to be more vulnerable to extinc- 
tion as evidence by recent extinctions (Roberts and Hawkins 
1999; Dulvy et al. 2003). These extinctions highlight the need 
to examine the genetic characteristics of island fish faunas 
to help determine why this group may be vulnerable. Pro- 
cesses that influence genetic resilience and promote species 
persistence times may be identified by examining spatial and 
temporal patterns of gene flow. 

Coral reef anemonefish (genus Amphiprion) represent a 
useful model system for understanding evolutionary histo- 
ries and population genetic structures of island fish faunas. 
Although the genus has a broad Indo-Pacific distribution, 
more than 25% of species are endemic to isolated islands or 
have peripheral populations at these remote locations (Fautin 
and Allen 1997). Although these island endemics are often 
closely related to more broadly distributed sister species, the 
historic colonization and speciation processes, and current 
levels of population differentiation are not well known. In a 
detailed phylogenetic study of 23 of the 28 anemonefish, 
Santini and Polacco (2006) suggested the group (Family: 
Pomacentridae, subfamily Amphiprioninae) originated some 
5-13 million years ago in the Indo-Pacific, with many of the 
endemics being of recent origin. A more detailed investiga- 
tion of their historic relationships with putative sister species, 
and current levels of gene flow among locations, is necessary 
to understand how endemic anemonefish originate and per- 
sist. This information will be useful for predicting how these 
species will persist in the future and aid management strate- 
gies aimed at conserving these iconic coral reef fish. 

Using phylogenetic and population genetic analyses, we 
reconstructed the evolutionary history of two Australian 
anemonefish: a small-range species, Amphiprion mccullochi 
(endemic to Middleton Reef [MR], Elizabeth Reef [ER], 
and Lord Howe Island [LHI]) and its more widespread sis- 
ter species, A. akindynos, found on the Great Barrier Reef, 
New Caledonia and the subtropical east coast of Australia. 
The endemic McCulloch's anemonefish is of particular con- 
servation concern because it has three characteristics known 
to elevate the risk of extinction — a very small geographic 
range (Coleman 1980; Fautin and Allen 1992; Hobbs et al. 
2009a), extreme habitat specialization (one species of host 
anemone — Entacmaea quadricolor; Fautin and Allen 1992) 
and relatively small local populations (Choat et al. 2006; 
Hobbs and Feary 2007). Entacmaea quadricolor is distributed 
from Micronesia and Melanesia to East Africa and the Red 
Sea and to Japan and Australia (Fautin and Allen 1992), but 
no further south than the Solitary islands (NSW) on the 
east coast of Australia. Amphiprion mccullochi is thought to 
have arisen by divergence from a more widespread most re- 
cent common ancestor (mrca) shared with A. akindynos, its 
sister species, on the eastern Australian coast (Santini and 
Polacco 2006). However, the mode of speciation and current 
levels of genetic connectivity among the three populations of 
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A. mccullochi are unknown. Amphiprion akindynos may be 
less of a conservation concern because it is more widely dis- 
tributed, inhabits six species of anemones (Fautin and Allen 
1997), and can be locally abundant (e.g., Richardson 1999). 

This study addressed the following specific questions: 
(1) What is the evolutionary relationship between the sis- 
ter species, A. mccullochi and A akindynos, based on mito- 
chondrial DNA (mtDNA)? (2) How many management units 
(MUs, sensu Moritz 1994) can be identified for A. mccullochi 
and A. akindynos 7 . (3) What is the contemporary relation- 
ship of connectivity between and within species, based on 
msat DNA? (4) What are the genetic diversities of these sis- 
ter species and do they suggest resilience or susceptibility to 
extinction and environmental change? Identifying MUs, in- 
cluding the direction of connectivity between these units, is 
essential information to ascertain best practice management 
and maximize biodiversity conservation of these important 
southernmost coral reefs. We discuss the implications of the 
underlying genetic structure to extinction risk and the con- 
servation of these remote populations. 

Materials and Methods 
Study system 

McCulloch's anemonefish (A. mccullochi) has the smallest 
geographic range of any of the 28 species of anemonefish 
and is endemic to three isolated, oceanic locations more than 
600 km off the east Australian coastline (ER, MR, and LHI). It 
is only found living in close association with its host anemone, 
E. quadricolor, and occurs at depths between 2 and 45 m 
(Fautin and Allen 1992). Its coloration of black body with 
a whitish snout, caudal peduncle, and caudal fin makes it 
easily recognizable. Its sister species, A. akindynos (Santini 
and Polacco 2006) is more widespread, ranging from the 
Great Barrier Reef (GBR) south to the Solitary Islands and 
extending out to New Caledonia (but not including ER and 
MR or LHI). It lives among its host anemones E. quadricolor, 
Heteractis aurora, H. crispa, H. magnifica, Stichodactyla 
haddoni, and S. mertensii (Fautin and Allen 1992), at depths 
between 1 and 25 m (Allen 1991). It has two white bars on its 
body and a color transition from a dark brown/orange body 
to a whitish caudal fin (Fautin and Allen 1992). 

Sampling locations and procedures 

Finclips from 60 A. mccullochi individuals were collected from 
two out of the three known populations at MR ( n = 30; Choat 
et al. 2006) and from LHI (n = 30; Hobbs et al. 2009a), and 
preserved in 70% alcohol. These two populations (MR and 
LHI) represent either end of the entire geographic range of 
this species, spanning 200 km. Bay et al. (2006) sampled 
A. akindynos at two GBR locations — a central, Lizard Island 
(LI) population ( n = 20) and a southern peripheral, One Tree 
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Island (OTI) population (n = 24), spanning 1200 km. The 
geographic ranges of the two sister species are separated by 
at least 600 km of deep open ocean habitat. 

Genetic techniques 

Standard DNA laboratory and analytical procedures were 
used for mtDNA sequencing (Bay et al. 2006, van Herwerden 
et al. 2009). Total genomic DNA was extracted from approx- 
imately 1 mm 3 of tissue using standard chelex-proteinase 
K digestion extraction procedures (Walsh et al. 1991). Each 
20-u.L polymerase chain reaction (PCR) amplification re- 
action contained 20 ng DNA, 0.2 mM DNTP, 1 unit 
of Bioline Biotaq Red DNA Polymerase, 0.5 /xM of 
each primer, and variable MgCl2. The noncoding re- 
gions of the mtDNA were selectively amplified by poly- 
merase chain reaction (PCR) using the following primers 
for the damselfish Acanthachromis polyacanthus: dLoop F 
(5'-CATATATGTRTTATCAACATTA-3') and CR-E H16498R 
(5'-CCTGAAGTAGGAACCAGATG-3') (Bay et al. 2006). 
Primers were tested and optimized using a Bio-Rad CI 000 
Thermal Cycler (Bio-Rad, Australia). Amplifications fol- 
lowed the same basic cycling protocol: 40 sec at 94° C, 
40 sec at primer-specific annealing temperatures (two dif- 
ferent touchdown profiles of five cycles at 55°C followed by 
30 cycles at 53°C or five cycles at 53°C followed by 30 cycles at 
5 1°C) and 40 sec at 72°C. The cycling profile was flanked by an 
initial 3-min denaturing step (94°C) and a 5-min terminal ex- 
tension phase (72°C). Furthermore, we tested and sequenced 
a 16S mitochondrial marker (Santini and Polacco 2006) on 
10 individuals per species (A. mccullochi and A. akindynos) . 
The resulting phylogeographic tree was poorly resolved and 
provided no resolution. Thus, only D-loop results are pre- 
sented in this study to identify haplotypic differences. We 
also genotyped both A. mccullochi and A. akindynos using 
17 microsatellite markers developed for A. mccullochi (van 
der Meer et al. 201 1). Six multiplex reactions were performed 
and each 10-uX PCR amplification reaction contained 10 ng 
DNA, 5 u.L (2x) Type-it Multiplex PCR Master Mix, 2 nL 
multiplex primer mix (at a concentration of 2 /xM). Ampli- 
fications followed the same basic cycling protocol: 30 sec at 
95°C, 90 sec at 58°C, and 40 sec at 72°C. The cycling profile 
was flanked by an initial 5-min denaturing step (94°C) and a 
30-min terminal extension phase (72°C). PCR products were 
purified by isopropanol precipitation for direct sequencing of 
D-loop (Macrogen, Korea), and ethanol ammonium acetate 
precipitation for genotyping microsatellite reactions (Genetic 
Analysis Facility, James Cook University, Townsville). 

mtDNA analysis 

Forward sequences were automatically aligned using the plu- 
gin CLUSTAL W in Geneious Pro 4.7 (Drummond et al. 
2009), conservatively trimmed to minimize the amount of 



missing data, and manually edited, inserting gaps where re- 
quired and checking for ambiguities. Differences between 
individual sequences were determined for the following 
characters: A, G, T, C, and IUB symbols (Nomenclature 
Committee 1985). Sequence data were obtained from Gen- 
Bank for the following five species that acted as outgroups: 
A. clarkii (DQ343928.1), A. clarkii (orange; DQ343929.1), 
A. chrysopterus (DQ343927.1), A. latezonatus (DQ343933.1), 
and A. leucokranos (DQ343934.1). All A. akindynos sequences 
(DQ250449.1 to DQ250492.1) from Bay et al. (2006) were 
analyzed together with A. mccullochi sequences based on the 
close genetic relationship between these two species (Santini 
and Polacco 2006). One GBR A. akindynos sequence from an 
unspecified location was also included (DQ343924.1). 

jModeltest (Posada 2008) identified an HKY + I model 
under Akaike Information Criterion with gamma = 0.271. 
Fifty- three of the 317 nucleotides sequenced for A. mccullochi 
were parsimony informative. The transition (ts): transver- 
sion (tv) substitution ratio was approximately 6:1. There 
were a total of 65 variable sites excluding five single base in- 
dels. The nucleotide composition was AT biased with 71.11% 
AT:28.89% GC (A. mccullochiDloop), 71.60% AT:28.40% GC 
(A. afomfynos Dloop),and71.36%AT:28.64%GC (combined 
A. mccullochi and A. akindynos Dloop), which is consistent 
with fish mitochondrial DNA (McMillan and Palumbi 1997). 

The data underwent four phylogenetic analyses (Fig. 1A): 
(1) Bayesian inference (MB) in MrBayes 3.1 (Huelsenbeck et 
al. 2001, 2002) with 10 million generations of Monte Carlo 
Markov chains (MCMC); (2) ten independent maximum 
likelihood (ML) analyses, followed by an 100 bootstrap repli- 
cate analysis using GARLI 0.951 (Zwickl 2006), from which 
a 50% majority rule ML consensus tree was constructed 
in PAUP* 4.10b (Swofford 2001); (3) maximum parsimony 
(MP) was performed in MEGA 4.0 (Tamura et al. 2007) 
with 1000 bootstrap replicates, from which a 50% majority 
rule consensus tree was constructed and; (4) Bayesian in- 
ference in BEAST VI. 6.1 (Drummond and Rambaut 2007) 
was tested using a strict clock model (estimated clock rate, 
uniform prior distribution, HKY+G site model for 5 million 
MCMC chains with sampling at every 5000 trees) and a re- 
laxed clock model (uncorrelated lognormal clock model, with 
the same parameters noted previously). Tracer VI. 5 identi- 
fied no significant difference between clock models based on 
the Bayes Factor (BF) evaluation of the models (Kass and 
Raftery 1995; Suchard et al. 2001), BF = 0. We therefore 
used the speciation (Yule process) to construct a Bayesian 
skyline plot using a strict clock model (as above, but with 
50 million MCMC chains sampled every 5000th tree). Max- 
imum clade credibility trees (MCCT) were constructed in 
TreeAnnotator VI. 6.1 after discarding the initial burn-in of 
10%. The MCCT was viewed separately in FigTree VI. 3.1 
(available at http://tree.bio.ed.ac.uk/software/figtree/). The 
best outgroup rooted ML tree from GARLI was selected to 
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Figure 1. (A) An outgroup rooted phylogram of D loop sequences from 50 Amphiprion mccullochi individuals from Middleton Reef and Lord Howe 
Island and 44 A. akindynos individuals from Lizard and One Tree Islands on the Great Barrier Reef. Amphiprion clarkii was specified as the outgroup 
species. This represents the best ML tree from 10 individual runs in GARLI. Asterisks identify individuals of A. akindynos that fall outside a clade. 
Numbers on branches indicate support for each clade, BEAST, MB, ML, and MP. (B) Haplotype minimum spanning tree (MST) with number of 
substitutions between haplotypes indicated on connectors. Different shading represents each of the four locations as shown on the key to the figure. 



reconstruct the evolutionary history with bootstrap values 
for each clade from all four analyses, if present (Fig. 1A). A 
minimum spanning tree (MST) of the mtDNA sequence data 
was computed using Arlequin ver. 3.5 (Excoffier et al. 2005) to 
identify shared haplotypes between locations sampled and/or 
species (Fig. IB). 

Bayesian skyline plots in BEAST VI. 6.1 were constructed 
to evaluate the presence of demographic stability and/or ex- 
pansions from coalescent analyses using strict clock model 
(parameters as above). We also evaluated population stasis 
using Fu's Fs parameter for population stasis (Fu 1997) and 
Tajima's D test for selective neutrality of mtDNA (Tajima 
1983), both of which accepted the hypothesis of a static 



population under an assumption of selective neutrality for 
A. mccullochi (D = 1.578, P = 0.96 and Fs = 2.80, P = 0.85); 
however, for A. akindynos selective neutrality of the mtDNA 
was accepted (D = 0.103, P = 0.61), while population 
stasis was rejected (Fs = -24.42, P < 0.0001), suggesting 
a spatial expansion by A. akindynos. Following 90,000 re- 
samplings of the data, the F -statistics (fixation indices, <t> sf , 
Od, <E> SC ) were determined using an analysis of molecular 
variance (AMOVA) (Excoffier et al. 1992) to detect popula- 
tion genetic partitioning between either regions (LI and OTI 
vs. MR and LHI) or genetically distinct lineages regardless 
of location or species (MUs were compared to each other; 
Table 1). Haplotype (h) and nucleotide diversities (%tt) of 
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Table 1. Analysis of molecular variance (AMOVA) analysis for the genetic data from Amphiprion mccullochi and A. akindynos structured into 
(A) mtDNA — partitioned by geographic region (GBR vs. offshore locations) and (B) mtDNA — partitioned by ESUs, and (C) msat DNA — partitioned by 
geographic region (GBR vs. offshore locations). 

F-statistics fixation indices 

Source of variation df SSD Variance component Percentage of variation (P-value) 

(A) Region 

Among groups 1 170.34 3.53 32.63 F a = 0.326 

(0.341 ± 0.014) 

Among populations within groups 2 10.10 -0.10 -0.93 F sc = -0.014 

(0.596 ±0.011) 

Within populations 90 665.51 7.39 68.31 F st = 0.317 

(<0.001 ± 0.000) 



(B) Clades 

Among groups 1 352.29 6.38 49.56 F a = 0.496 

(0.007 ± 0.003) 

Among populations within groups 10 367.64 5.13 39.88 F sc = 0.791 

(<0.001 ± 0.000) 

Within populations 80 108.74 1.36 10.56 F st = 0.894 

(<0.001 ± 0.000) 

(C) Microsatellite 

Among groups 1 30.54 0.25 5.01 F a = 0.05 

(0.337 ± 0.002) 

Among populations within groups 2 13.684 0.04 0.85 F sc = 0.009 

(0.006 ± 0.000) 

Within populations 192 917.64 4.78 94.14 F st = 0.05 

(<0.001 ± 0.000) 



the data were interpreted based on Grant and Bowen (1998). 
Population pairwise F st comparisons (measured in Arlequin 
ver. 3.5) initially identified no differences between locations 
for either species, which informed further inter- and intraspe- 
cific analyses between clades (Table 2). Isolation by distance 
(IBD) between reefs was tested using a Mantel test in IBD 
verl.4 with 10,000 permutations (Bohonak 2002). 

Microsatellite analyses 

GENEPOP 4.0 (Rousset 2008) was used to perform exact tests 
of departures from Hardy-Weinberg equilibrium (HWE) for 
each locus per sampled location (i.e. 17 loci x four sam- 
pled locations) and to test for linkage disequilibrium (LD) 
between the 17 loci within each of the two study species 
(i.e., 17 x 17 (-17) = 272 tests for each species), using the 
Markov chain algorithm. If departure from HWE was ob- 
served, the program MICRO-CHECKER 2.2.3 (van Ooster- 
hout et al. 2004) was used to detect the presence of null alleles, 
large allele dropout, and other scoring errors. We conducted 
20 batches with 5000 iterations per batch. A false discov- 
ery rate (FDR; Benjamini and Hochberg 1995) correction 
was applied to all HWE and LD results, using the program 
QVALUE (Storey 2002). Significant single-locus departures 
from HWE were detected in eight of sixty-eight tests at the 
population level before FDR correction and five afterwards 



(Ami, Amll, Aml4, Aml7, Aml9). Null alleles were iden- 
tified in LI (Ami 1, Aml7, Aml9), OTI (Ami, Am7, Am21), 
MR (Amll,Aml4,Aml7),andLHI (Amll, Aml7, Aml9) as 
indicated in MICROCHECKER. Loci that were not in HWE 
and had null alleles (i.e., Ami, Amll, Aml4, Aml7, Aml9) 
were not used in subsequent analysis (ARLEQUIN, STRUC- 
TURE, NEWHYBRIDS, and MIGRATE-n). Of 544 locus x 
locus exact tests of LD, (272 per species), only 13 were signif- 
icant before FDR and none after FDR correction (Benjamini 
and Hochberg 1995), indicating that loci are independently 
assorted. 

Following 90,000 resamplings of the data, P-statistics (fixa- 
tion indices, O st , Oct) ^sc) were determined using an AMOVA 
(Excoffier et al. 1992) to detect population genetic parti- 
tioning between regions (LI and OTI vs. MR and LHI; 
Table 1). Microsatellite gd of differentiated populations was 
determined in ARLEQUIN 3.5 (Excoffier et al. 2005), using 
90,000 permutations (Table 2). IBD was tested using a Man- 
tel Test in IBD verl.4 with 10,000 permutations (Bohonak 
2002), as noted for mtDNA analyses. 

We used STRUCTURE V2.3 (Pritchard et al. 2000; Hubisz 
et al. 2009) to identify groups of randomly mating individu- 
als based on microsatellite allele frequencies. Data were first 
tested using the "admixture" and "no admixture" models, 
including information about the sampling location as a prior 
and correlated allele frequencies. Two repetitions were run 
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Table 2. Genetic diversity estimates for Amphiprion mccullochi and A. akindynos. Sample size (n), number of haplotypes (n h ), haplotypes diversity ± 
SE (h), nucleotide diversity ± SE (%jt), and genetic diversity (gd). Genetic diversity estimates for one A. akindynos sample of unknown location on the 
GBR was omitted. 



Site 


n 


ih 


h 


%7T 


gd 


A. mccullochi 












Middleton Reef 


22 


14 


0.952 ± 0.026 


4.90 ± 2.53 


0.691 ±0.360 


Lord Howe Island 


26 


1 1 


0.846 ± 0.054 


5.1 1 ± 2.62 


0.670 ± 0.351 


Total 


48 


21 


0.905 ± 0.027 


5.00 ± 2.50 


0.690 ± 0.370 


Clade 1.1 


22 


9 


0.701 ± 0.103 


1.76 ±0.98 




Clade 1.2 


6 


3 


0.600 ± 0.215 


0.1 1 ± 0.14 




Clade 2.4 


1 1 


4 


0.600 ±0.154 


0.22 ± 0.20 




Clade 2.5 


3 


3 


1.000 ±0.272 


0.63 ± 0.60 




Clade 2.6 


8 


4 


0.643 ±0.184 


0.24 ± 0.22 




A. akindynos 












Lizard Island 




1 A 
I 4 


u.yi>z ± u.uzd 


A OA -1- T QZt 

4.yu ± Z.dd 


n 71 q -i- n Q"7 1 
U. / I o ± U.J / I 


One Tree Island 


24 


11 


0.846 ±0.054 


5.11 ±2.62 


0.693 ±0.360 


Total 


44 


34 


0.976 ± 0.013 


4.03 ± 2.06 


0.706 ± 0.484 


Clade 1.3 


5 


5 


1.000 ±0.127 


1 .46 ± 1 .00 




Clade 1.4 


8 


4 


0.643 ±0.184 


0.40 ±0.32 




Clade 2.1 


2 


2 


1.000 ±0.500 


0.32 ±0.45 




Clade 2.2 


3 


3 


1.000 ±0.272 


0.63 ±0.60 




Clade 2.3 


8 


6 


0.893 ±0.1 11 


0.76 ±0.52 




Clade 2.7 


5 


4 


0.900 ±0.161 


1.20 ±0.85 




Clade 2.8 


1 1 


7 


0.81 8 ± 0.1 19 


0.56 ±0.40 




A. mccullochi and A, akindynos 












Total 


93 


55 


0.969 ± 0.009 


5.67 ±2.81 


0.691 ±0.370 



for each possible K value, informed by the number of parti- 
tions identified by mtDNA analyses, so K ranged between 1 
and 12. Short 20,000 MCMC iterations were performed per 
analysis and a 1000 tree burn-in was applied. To determine 
the "best value" for K, we followed the method suggested 
by Pritchard et al. (2000), which involved comparing mean 
log likelihoods penalized by one-half of their variances (see 
Hubisz et al. 2009). We also chose an "admixture" model in- 
cluding information about the sampling location as a prior 
and correlated allele frequencies, with two repetitions of K 
(K = 2), for 1 million MCMC iterations and a 10,000 tree 
burn-in, as this was the best K value identified in the initial 
analyses. 

Given contrasting results found with mtDNA and msat 
DNA, we specifically tested for interbreeding using two 
programs: (1) NEWHYBRIDS (Anderson and Thompson 
2002), which implements a Bayesian method aimed at de- 
tecting the presence of hybrids from a sample of individ- 
uals of mixed origin. We used a MCMC procedure with 
a 150,000 burn-in and 150,000 steps, and (2) MIGRATE- 
n 2.4.3 (http://popgen.sc.fsu.edu/Migrate-n.html; see Beerli 
and Felsenstein 2001; Beerli 2004) to estimate long-term gene 
flow between MUs and short-term gene flow (nuclear DNA) 
between locations. We set the migration rate parameter for 
mtDNA (Theta and M to a maximum of 0. 1 and 5000, respec- 



tively) and msatDNA (both Theta and M to a maximum of 
100). We conducted 10 replicates of a Bayesian analysis with 
one long-chain sampling every 100 trees of 100,000 sampled 
and a 20,000 iteration burn-in for mtDNA and a Bayesian 
analysis with one long-chain sampling every 100 trees of 
5000 sampled and a 1000 iteration burn-in for msatDNA. All 
parameters converged and fell within the 90% CI, yielding 
values for 0 and m for each locus per population. Finally, a 
discriminant analysis of principal components (DAPC) was 
also used to determine contemporary gene flow between pop- 
ulations as per Jombart et al. (2010). This required the use of 
the program R 2.12 (http://www.Rproject.org). We retained 
75 principal components (PCs) comprising 95% of the total 
genetic information as predictors for DA. 

Results 

Evolutionary relationship 

The best outgroup rooted ML phylogram and the MST identi- 
fied two strongly supported evolutionary groups and this was 
supported by all additional phylogenetic analyses (Fig. 1 A). It 
appeared that either one evolutionary group emerged from 
the other or that there was a splitting to form two sister 
groups, depending on the analysis. The two evolutionary 
groups did not represent the two species per se, and lacked 
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reciprocal monophyly, because individuals of both species 
were detected in both evolutionary groups. Despite the lack 
of species-specific partitions, there were no shared haplo- 
types between species. This suggests a complex evolution- 
ary history between these recognized species or incomplete 
lineage sorting. Sixteen synapomorphic substitutions char- 
acterized the split between these two phylogenetic groups ( 1 
and 2). Group 1, which emerges from Group 2, contains 56% 
of all sampled A. mccullochi, whereas the apparently older 
Group 2 contains 70% of all sampled A. akindynos repre- 
sented. 

Management units 

The two evolutionary groups contained a total of seven and 
possibly as many as 12 MUs, of which three to four were 
in Group 1 (MU1-4) and four to eight (MU5-12) were in 
Group 2. Haplotype sharing was observed between loca- 
tions within species: MR — LHI in all five A. mccullochi MUs; 
LI — OTI in five of the seven A. akindynos MUs. When sam- 
ples were structured by MU irrespective of sampling loca- 
tion (MU1-12) during population genetic analyses, pair- 
wise F st values indicating significant genetic differentiation 
were obtained between all but two MU pairs (pairwise 
F st = 0.455-0.983, P < 0.00001 to 0.04800). Exceptions 
were MU5-6 (pairwise f st = 0.634, P = 0.10) and MU5-9 
(F st = 0.860, P = 0.10). Lack of significance among these 
MU-pairs, despite very large pairwise F st values may result 
from smaller sample sizes in these three clades than in the re- 
maining nine clades (see Fig. 1A). AM OVA results confirmed 
genetic partitioning between MUs, O st = 0.894 (P < 0.001); 
however, this only explained 10.56% of the variation, while 
49.56% and 39.88% of the variation occurred among groups 
(1 vs. 2) and among populations within groups, respectively 
(<J> ct = 0.496, P < 0.001; 4> sc = 0.791, P < 0.001; Table 1). 

We further examined levels of gene flow between loca- 
tions and/or species using pairwise F s ^ values for a number 
of interspecies and mixed-species comparisons, consistent 
with the phylogenetic signals described above. There was no 
mtDNA differentiation between sample locations (LHI vs. 
MR or OTI vs. LI) within either species (A. mccullochi or 
A. akindynos, F st = -0.005 and -0,025, P = 0.414 and 0.865, 
respectively). Given this and disregarding species identity, to 
compare fish from offshore to fish from continental shelf 
locations (OTI or LI vs. LHI or MR), all combinations of 
comparisons showed significant pairwise mtDNA F st values 
(F st = 0.293-0.349, P < 0.00001). The more statistically rig- 
orous AMOVA of samples structured by region (LHI and 
MR vs. OTI and LI) confirmed mtDNA genetic partitioning 
with more than two-thirds (68.31%) of the genetic variation 
within locations, 4> st = 0.317 (P < 0.0001), and less than 
one-third of the variation detected among regions (offshore 
vs. GBR continental locations), which was not significant, 



<J) ct = 0.327 (P = 0.341). Further, none of the variation oc- 
curred among locations within regions (i.e., LHI vs. MR; OTI 
vs. LI), <I> SC = -0.014 (P = 0.596). There was no IBD based 
on mtDNA from all locations sampled, using a Mantel test of 
pairwise geographic (km) and genetic (F st ) distance between 
locations (z = 2079.08, R 2 = 0.31, P = 0.337). 

Both Fu's Fs parameter for population stasis (Fu 1997) and 
Tajima's D test for selective neutrality of mtDNA (Tajima 
1983) accepted the hypothesis of a static population un- 
der an assumption of selective neutrality for A. mccullochi 
(D = 1.578, P = 0.96 and Fs = 2.80, P = 0.85). In contrast, 
population stasis was rejected for A. akindynos (D = -24.42, 
P < 0.0001; Fs = -11.757, P = 0.02), suggesting that this 
species has undergone spatial expansion. This was confirmed 
by both Bayesian and mismatch distribution analyses (data 
not shown). When mtDNA data of specimens from all loca- 
tions, regardless of species partition, were considered collec- 
tively (MR vs. LHI vs. OTI vs. LI), total mtDNA suggested 
selective neutrality (D = 0.103-0.636, P = 0.61-0.80) and 
appeared to be either evolutionarily stable or to represent 
admixture of previously differentiated lineages. 

Contemporary connectivity 

Msat pairwise F st values were largely consistent with 
mtDNA results: There was no genetic differentiation between 
A. mccullochi populations, LHI-MR (F st = 0.005, P = 0.081), 
but significant differentiation existed between A. akindynos 
populations, LI and OTI (F st = 0.016, P = 0.036). Impor- 
tantly, the more statistically rigorous AMOVA of samples 
structured by region (LHI and MR vs. OTI and LI) identi- 
fied msat genetic partitioning within locations, accounting 
for 94.14% of the genetic variation; <J> st = 0.06 (P < 0.001). 
Only 5.01% of the variation occurred among regions (off- 
shore vs. GBR continental locations), but this was not signif- 
icant, O ct = 0.05 (P= 0.337), and almostnone (0.85%) ofthe 
variation, albeit significant, occurred among locations within 
regions (i.e., LHI vs. MR, OTI vs. LI), O sc = 0.009 (P= 0.006; 
Table 1). In contrast to mtDNA, there was strong evidence of 
IBD based on msat DNA using a Mantel test of log pairwise ge- 
ographic (km) and log genetic (F st ) distance between all sam- 
pled locations, regardless of species (z = -23.11, R 2 = 0.89, 
P = 0.048). 

NEWHYBRIDS and STRUCTURE analyses indicated in- 
traspecific but not interspecific gene flow between species 
and the likelihood of the marginal posterior probability dis- 
tribution was greatest when K = 2 (Fig. 2A and B). DAPC 
was consistent with this and A. mccullochi populations (MR 
and LHI) appeared to have overlapping genotypic profiles, as 
per mtDNA, while A. akindynos populations (OTI and LI on 
the GBR) were genetically distinct from each other (Fig. 2C), 
unlike mtDNA results. Using the four sampled populations 
as a priori population criteria, DAPC assigned 96.94% of 
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Figure 2. Separation of Amphiprion mccullochi and A. akindynos based on various analyses of msat loci: (A) NEWHYBRID analysis showing pure 
(gray) and F1 hybrid (black) status, (B) STRUCTURE cluster analysis, and (C) Scatterplots of the discriminant analysis of principal components (DAPC) of 
msat data for two A. mccullochi and two A. akindynos populations using geographic sample site as priors for genetic clusters. Populations are named 
and individual genotypes appear as dots surrounded by 95% inertia ellipses. Eigenvalues show the amount of genetic information contained in each 
successive principal component with X and Y axes constituting the first two principal components, respectively. 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1599 



Historic Hybridization between Anemonefish 

all individuals to the population where they were sampled 
(assignment per population, MR = 80%, LHI = 80%, 
LI = 100%, and OTI = 100%). MIGRATE-n indicated high 
levels of historical gene flow relative to contemporary gene 
flow. Levels of historical gene flow between lineages (i.e., 
Group 1 vs. Group 2; 4 Nm values ranged from 28 to 62) 
were an order of magnitude lower compared to within lin- 
eage geneflow: Group 1 (MU1. 1-1.4; 4 Nm values ranged 
from 980 to 2326) and Group 2 (MU2. 1-2.8; 4 Nm values 
ranged from 241 to 1247). Contemporary gene flow was a 
few orders of magnitude less than historical gene flow, with 
4 Nm values ranging from 1 to 6. 

Genetic diversity 

Genetic diversity was high for both species, evident from 
high haplotype (h), nucleotide (%n), and genotypic (gd) 
diversities for each sampled location (h = 0.846-0.982, 
%7i = 3.97-5.11, gd = 0.670-0.718, for each species). 
Total haplotype, nucleotide, and genotypic diversity were 
also high when species were combined (h = 0.905-0.976, 
%n = 4.03-5.67, gd = 0.690-0.706; Table 2). When samples 
were grouped according to mtDNA evolutionary groups, ir- 
respective of species or location, haplotype diversity was still 
high (h = 0.600-1.000), but nucleotide diversity was at least 
one-third (%ir = 0.11-1.76) that of the aforementioned val- 
ues for each species (Table 2). 

Discussion 

Phylogenetic analysis of mtDNA revealed a complex evolu- 
tionary history where A. mccullochi and the more widespread 
A. akindynos were not partitioned into monophyletic lineages 
as anticipated, but were mixed into two paraphyletic lineages 
that did not correspond to either species or geographic lo- 
cation per se. These analyses also revealed that the species 
did not share any haplotypes with each other and identified 
up to 12 species-specific MUs distributed among the four 
geographic locations sampled for the two species. Population 
genetic analyses based on mtDNA confirmed the absence of 
genetic partitioning by location for both species examined. 
Importantly, despite this complex evolutionary history, msat 
DNA analyses revealed no contemporary gene flow existed 
between species, and geneflow between populations was lim- 
ited. Populations connected by such low levels of gene flow 
are effectively demographically independent and significant 
to conservation as these MUs are fundamental to effective 
demographic management and are the logical unit for de- 
mographic study and population monitoring (Moritz 1994). 
Finally, the high genetic diversity of both species may result 
in greater adaptive capacity to cope with future environmen- 
tal change, insofar as greater genetic diversity provides more 
"raw material" for selection to act on. 
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Evolutionary relationship 

There are three possible interpretations of the combined ge- 
netic results. First, the two species represent different color 
variants of a single species rather than different species per se. 
However, this contradicts existing taxonomic classification 
(Santini and Polacco 2006) and requires further ecological 
and experimental data to explore this. The second possibility 
is that incomplete lineage sorting may be responsible for the 
apparent lack of reciprocal mtDNA monophyly. However, 
as mtDNA lineages suffer incomplete lineage sorting for a 
much shorter period of time (25%) than do nuclear DNA 
lineages, our msat DNA results do not support this scenario, 
because msat DNA did partition the species, which leads to 
the third possible interpretation — that there has been histor- 
ical bi-directional hybridization. 

We argue that the high level of mtDNA genetic partition- 
ing observed — while not identifying species-specific lineages 
or location-specific MUs — is most consistent with an evo- 
lutionary history of at least two reticulate events between 
A. akindynos and A. mccullochi. Reticulate events have been 
documented for numerous other reef fish species (e.g., 
McMillan et al. 1999; van Herwerden and Doherty 2006; 
van Herwerden et al. 2006; Yaakub et al. 2006; Marie et al. 
2007), including secondary contact between differentiated 
lineages of A. akindynos on range edges as documented by 
Bay and Caley (20 1 1 ) . Hybridization and reticulate evolution 
is also common in other coral reef organisms (e.g., corals: 
Willis et al. 2006 ) . Several lines of evidence support a scenario 
of historical hybridization: (1) two paraphyletic mtDNA lin- 
eages exist, each comprising of both A. mccullochi and A. 
akindynos, where each lineage may represent a mtDNA lin- 
eage of one of the species prior to historic hybridization; 
(2) there are no shared mtDNA haplotypes between present 
day populations of the two species; demographic geneflow is 
very limited between the species, as measured by msat DNA 
and Migrate-n; (3) there is population genetic partitioning 
between the two species, evident from both mtDNA and 
msatDNA, even if treated as populations rather than species 
and; (4) msat genotypes of almost 97% of individuals were 
assigned back to the population they were collected from and 
there was strong and total partitioning between species fol- 
lowing DAPC analysis of msat DNA. Together, this suggests 
that geneflow between species occurred in the evolutionary 
past but is either no longer occurring, or is happening at a 
level not detectable in this study. 

Historical gene flow between species via hybridization is 
increasingly being documented among coral reef fish (e.g., 
Yaakub et al. 2006; Hobbs et al. 2009b, in press) and has 
been reported in a number of anemonefish species, including 
A. akindynos (Fautin and Allen 1997; Sea Read 2009) . Timm et 
al. (2008) reports hybridization between Amphiprion species, 
as a possible explanation for sequence sharing, especially 
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since, within the genus Amphiprion, several species have sim- 
ilar coloration and overlapping variation at otherwise diag- 
nostic morphological characters. Coloration is very pertinent 
in the present case, as juvenile A. akindynos and A. mccullochi 
are almost morphologically indistinguishable (Richardson 
1998). Settlement to an anemone occurs during the juvenile 
phase and given that both species use the same host anemone 
(E. quadricolor; Fautin and Allen 1992) and have similar col- 
oration, this could lead to the formation of heterospecific 
social groups, and possible interbreeding. Furthermore, al- 
though the two species occur in allopatry, two vagrant in- 
dividuals of A. akindynos have been recorded at LHI (Crean 
et al. 2009). The arrival of vagrants into the distributional 
range of an allopatric sister species is thought to promote 
interbreeding between species due to the low availability of 
conspecific partners for the vagrants (Hobbs et al. 2009b). 
The East Australian Current flows through the ranges of the 
two study species and the strength and direction is influ- 
enced by climate (Philander et al. 1990; Middleton et al. 
2009). Historical changes in this current may have facilitated 
the arrival of vagrants, which may have resulted in contact 
and hybridization between the two species. Taken together, 
all these factors suggest hybridization is a likely scenario in 
the evolutionary history of these two species. 

Management units 

Two groups (evolutionary units) and at least seven MUs were 
identified in this study, with at least two A. akindynos MUs 
and five A. mccullochi MUs. While each MU represents one 
of the study species exclusively, most are not partitioned 
by geography. This complex underlying phylogenetic struc- 
ture may have occurred because anemonefish have both self- 
recruitment at demographic scales (Jones et al. 2005) and in- 
terpopulation connectivity at longer evolutionary timescales 
(Timm et al. 2008). 

Contemporary connectivity 

Microsatellite loci detected very limited geneflow between 
the sister species A. akindynos and A. mccullochi, which 
strengthens our suggestion of evolutionary hybridization as 
current day populations are not mixing. Likewise, LI and OTI 
A. akindynos populations were genetically distinct which is 
consistent with Bay and Caley (201 1), indicating that the dis- 
tance between these two populations may be too great for 
geneflow to occur. In contrast, there was a lack of genetic 
partitioning for A. mccullochi populations. The difference 
between species in the level of gene flow between popula- 
tions probably represents different geographic distances be- 
tween the sample populations (McCulloch populations were 
160 km apart and the akindynos sample populations were 
1200 km apart). This gene flow between populations is 
promising for A. mccullochi since it facilitates recolonization 
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if one population was to go locally extinct. However, only two 
locations in each species were sampled and further samples 
of A. akindynos from the Great Barrier Reef, southern reefs 
(e.g., Solitary Islands), and New Caledonia, as well as samples 
of A. mccullochi from Elizabeth Reef will be needed to fully 
quantify gene flow between all populations of these species. 

Genetic diversity 

In both species, mtDNA data showed high genetic variabil- 
ity, h almost double and jt close to an order of magnitude 
greater than the cut-off defined by Grant and Bowen (1998; 
both h and jt >0.5). Similarly, high levels of genetic diversity 
for msatDNA were also found in both species. This suggests 
that populations of both species are either large and stable 
with long evolutionary histories, or that there has been sec- 
ondary contact between differentiated lineages (Grant and 
Bowen 1998). Given that A. mccullochi does not have a large 
population (Choat et al. 2006; Hobbs and Feary 2007), sec- 
ondary contact between differentiated lineages is more likely 
and is consistent with reticulate evolution between these sis- 
ter species. Given that mtDNA diversity tracks with nuclear 
genetic diversity in many marine species (reviewed by Jo- 
hannesson and Andre 2006), such high genetic diversity in 
A. mccullochi and A. akindynos is encouraging as it suggests 
that both species may have a greater adaptive capacity to deal 
with environmental change than if they had low genetic di- 
versity. However, a cautious approach is still warranted given 
that quantitative trait loci under selection, can have no ge- 
netic diversity in peripheral populations despite high neutral 
genetic diversity (Kellermann et al. 2009). 

Conserving genetic diversity is considered a priority by 
the IUCN (McNeely et al. 1990), and even more so for 
A. mccullochi given its vulnerability to extinction due to other 
traits (low abundance, small geographic range, and ecologi- 
cal specialization) . The high genetic diversity of A. mccullochi 
appears to have arisen through complex reticulate evolution 
involving historical (but not contemporary) hybridization. 
High genetic diversity and population connectivity may re- 
duce extinction risk in A. mccullochi, but identifying and 
protecting populations implicated in generating the com- 
plex reticulate structure among this species should be a con- 
servation priority. Conserving A. mccullochi would best be 
achieved by protecting each of its MUs and by minimizing 
threats to population size, such as habitat loss (e.g., anemone 
bleaching) and collection for the aquarium trade. 

Acknowledgments 

We thank the Australian Department of the Environment and 
Water Resources, and Envirofund Australia for funding. We 
also thank Lord Howe Island Marine Park Authority, Ian Kerr, 
Sallyann Gudge, and the captain and crew of the Capricorn 
Star for excellent logistical support. 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1601 



Historic Hybridization between Anemonefish 



M. H. van der Meer era/. 



References 

Allen, G. R. 1991. Damselfishes of the world. Mergus Publishers, 

Melle, Germany, p. 271. 
Allen, G. R. 2008. Conservation hotspots of biodiversity and 

endemism for Indo-Pacific coral reef fishes. Aquat. Conserv. 

18:541-556. 

Anderson, E. C., and E. A. Thompson. 2002. A model-based 

method for identifying species hybrids using multilocus 

genetic data. Genetics 160:1217-1229. 
Avise, J. C. 1992. Molecular population structure and the 

biogeographic history of a regional fauna: a case history with 

lessons for conservation biology. Oikos 63:62-76. 
Bay, L. K., and M. J. Caley. 201 1. Greater genetic diversity in 

spatially restricted coral reef fishes suggests secondary contact 

among differentiated lineages. Diversity 3:483-502. DOI: 

10.3390/d3030483. 
Bay, L. K., R. H. Crozier, and M. J. Caley. 2006. The relationship 

between population genetic structure and pelagic larval 

duration in coral reef fishes on the Great Barrier Reef. Mar. 

Biol. 149:1247-1256. 
Benjamini, Y., and Y. Hochberg. 1995. Controlling the false 

discovery rate, a practical and powerful approach to multiple 

testing. J. R. Stat. Soc. 57:289-300. 
Beerli, P. 2004. Effect of unsampled populations on the 

estimation of population sizes and migration rates between 

sampled populations. Mol. Ecol. 13:827-836. 
Beerli, P., and J. Felsenstein. 2001. Maximum likelihood 

estimation of a migration matrix and effective population sizes 

in n subpopulations by using a coalescent approach. Proc. 

Natl. Acad. Sci. 98:4563-4568. 
Bohonak, A. L 2002. IBD (Isolation By Distance): a program for 

analyses of isolation by distance. J. Hered. 93:153- 

154. 

Ceballos, G, and L H. Brown. 1995. Global patters of 

mammalian diversity, endemism and endangerment. Conserv. 

Biol. 9:559-568. 
Choat, J. H., L. van Herwerden, W. D. Robbins, J. P. Hobbs, and 

A. M. Ayling. 2006. A report on the ecological surveys 

conducted at Middleton and Elizabeth Reefs, February 2006. 

Report to the Australian Government Department of 

Environment and Heritage, Canberra, p. 65. 
Coleman, N. 1980. Australia sea fishes South of 30°S. Doubleday 

Australia Pty Ltd, Lane Cove, NSW, p. 309. 
Crean, A. J., S. E. Swearer, and H. M. Patterson. 2009. Larval 

supply is a good predictor of recruitment in endemic but not 

non-endemic fish populations at a high latitude coral reef. 

Coral Reefs 29:137-143. 
Drummond, A. J., and A. Rambaut 2007. Bayesian evolution 

analysis by sampling trees. BMC Evol. Biol. 7:214. DOI: 

10.1186/1471-2148-7-214 
Drummond, A. J., B. Ashton, M. Cheung, J. Heled, M. Kearse, R. 

Moir, S. Stones-Havas, T. Thierer, and A. Wilson. 2009. 

Geneious v4.7. Available at: http://www.geneious.com/ 



Dulvy, N. K., Y. Sadovy, and J. D. Reynolds. 2003. Extinction 

vulnerability in marine populations. Fish Fish. 4:25-64. 
Excoffier, L., P. E. Smouse, and ]. M. Quattro. 1992. Analysis of 

molecular variance inferred from metric distances among 

DNA haplotypes: application to human mitochondrial DNA 

restriction data. Genetics 131:479-491. 
Excoffier, L., G. Laval, and S. Schneider. 2005. Arlequin (version 

3.0), an integrated software package for population genetics 

data analysis. Evol. Bioinform. Online 1:47-50. 
Faith, D. P. 1992. Conservation evaluation and phylogenetic 

diversity. Conservation 61:1-10. 
Fautin, D. G., and G. R. Allen. 1992. Anemone fishes and their 

host sea Anemones: a guide for aquarists and divers. Western 

Australian Museum, Perth. 
Fautin, D. G., and G. R. Allen. 1997. Anemone fishes and their 

host sea Anemones: a guide for aquarists and divers. 2nd ed. 

Western Australian Museum, Perth. 
Frankham, R. 1997. Do island populations have less genetic 

variation than mainland populations. Heredity 78:311- 

327. 

Frankham, R. 1998. Inbreeding and extinction, Island 

populations. Conserv. Biol. 12:665-675. 
Fu, Y.-X. 1997. Statistical tests of neutrality of mutations against 

population growth, hitchhiking and background selection. 

Genetics 147:915-925. 
Gillespie, R. G., and G. K. Roderick. 2002. Arthropods on islands, 

colonization, speciation and conservation. Annu. Rev. 

Entomol. 47:595-632. 
Gillespie, R. G, E. M. Claridge, and G. K. Roderick. 2008. 

Biodiversity dynamics in isolates island communities: 

interaction between natural and human-mediated processes. 

Mol. Ecol. 17:45-57. 
Grant, W. S., and B. M. Bowen. 1998. Shallow population 

histories in deep evolutionary lineages of marine fishes: 

insights from the sardines and anchovies and lessons for 

conservation. ]. Hered. 89:415-426. 
Hobbs, ]. P., and D. A. Feary. 2007. Monitoring the ecological 

status of Elizabeth and Middleton Reefs, February 2007. 

Report to The Department of Water and Resources, Canberra, 

p. 37. 

Hobbs, ]. P., }. Neilson, and ]. ]. Gilligan. 2009a. Distribution, 
abundance, habitat association and extinction risk of marine 
fishes endemic to the Lord Howe Island region, November 
2009. Report to Lord Howe Island Marine Park, NSW Marine 
Parks Authority, p. 37. 

Hobbs, ]. P., A. J. Frisch, G. R. Allen, and L. Van Herwerden. 
2009b. Marine hybrid hotspot at Indo-Pacific biogeographic 
border. Biol. Lett. 5:258-261. 

Hubisz, M. J., D. Falush, M. Stephens, and J. K. Pritchard 2009. 
Inferring weak population structure with the assistance of 
sample group information. Mol. Ecol. Resour. 9:1322-1332. 

Huelsenbeck, ]. P., F. Ronquist, R. Nielsen, and J. P. Bollback. 
2001. Bayesian inference of phylogeny and its impact on 
evolutionary biology. Science 294:2310-2314. 



1602 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



M. H. van der Meer eta/. 



Historic Hybridization between Anemonefish 



Huelsenbeck, J. P., B. Larget, R. E. Miller, and F. Ronquist. 2002. 

Potential applications and pitfalls of Bayesian inference of 

phylogeny. Syst. Biol. 51:673-688. 
Johannesson, K., and C. Andre. 2006. Life on the margin: genetic 

isolation and diversity loss in a peripheral marine ecosystem, 

the Baltic Sea. Mol. Ecol. 15:2013-2029. 
lombart, T., S. Devillard, and F. Balloux. 2010. Discriminant 

analysis of principal components, a new method for the 

analysis of genetically structured populations. BMC Genet. 

11:94. 

Jones, G. P., M. J. Caley, and P. L. Munday. 2002. Rarity in coral 
reef fish communities. Pp. 81-102 in P. F. Sale, ed. Coral reef 
fishes. Dynamics and diversity in a complex ecosystem. 
Academic Press, San Diego. 

Jones, G. P., S. Planes, and S. R. Thorrold. 2005. Coral reef larvae 
settle close to home. Curr. Biol. 15:1314-1318. 

Kass, R. E„ and A. E. Raftery. 1995. Bayes factors. J. Am. Stat. 
Assoc. 773:795. 

Kellermann, V., B. van Heerwaarden, C. M. Sgro, and A. A. 
Hoffmann. 2009. Fundamental evolutionary limits in 
ecological traits drive Drosophila species distributions. Science 
325:1244-1246. 

Lessios, H. A., B. D. Kessing, and J. S. Pearse. 2001. Population 
structure and speciation in tropical seas, global 
phylogeography of the sea urchin Diadema. Evolution 
55:955-975. 

Marie, A. D., L. van Herwerden, J. H. Choat, and J. P. Hobbs. 

2007. Hybridisation of reef fish at the Indo-Pacific 

biogeographic barrier, a case study. Coral Reefs 26:841-850. 
McMillan, W. O., and S. R. Palumbi. 1997. Rapid rate of 

control-region evolution in Pacific butterflyfishes 

(Chaetodontidae). J. Mol. Evol. 45:473-484. 
McMillan, W. C, L. A. Weigt, and S. R. Palumbi. 1999. Color 

pattern evolution, assortative mating, and genetic 

differentiation in brightly colored butterflyfishes 

(Chaetodontidae). Evolution 53:247-260. 
McNeely, J. A., K. R. Miller, W. V. Reid, R. A. Mittermeier, and T. 

B. Werner. 1990. Conserving the world's biological diversity. 

World Conservation Union, World Resources Institute, 

Conservation International, World Wildlife Fund-US, and the 

World Bank, Washington, DC. 
Middleton, J., M. Baird, D. Cox, P. Estrade, J. Everett, A. Fowler, 

K. Hill, R. McCabe, T. Mullany, G. Nippard, et al. 2009. 

Topographic effects of the East Australia Current outflow at 

Lord Howe Island, Middleton Reef and Taupo Seamount. RV 

Southern Surveyor voyage SS06/2009 Voyage Plan and Voyage 

Summary. CSIRO Division of Marine and Atmospheric 

Research, Hobart, Tasmania. 
Moritz, C. 1994. Defining 'evolutionarily significant units' for 

conservation. Trends Ecol. Evol. 9:373-375. 
Moritz, C. 2002. Strategies to protect biological diversity and the 

evolutionary processes that sustain it. Syst. Biol. 51:238-254. 
Muss, A., D. R. Robertson, C. A. Stepien, P. Wirtz, and B. W. 

Bowen 2001. Phylogeography of Ophioblennius, the role of 



ocean currents and geography in reef fish evolution. Evolution 
55:561-572. 

Philander, S. G., J. R. Holton, and R. Dmowska. 1990. El Nino, La 

Nina, and the Southern Oscillation, Vol. 46 (International 

Geophysics). Academic Press, Waltham, 

Massachusetts. 
Planes, S., and C. Fauvelot 2002. Isolation by distance and 

vicariance drive genetic structure of a coral reef fish in the 

Pacific Ocean. Evolution 56:378-399. 
Posada, D. 2008. jModelTest, phylogenetic model averaging. Mol. 

Biol. Evol. 25:1253-1256. 
Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference of 

population structure using multilocus genotype data. Genetics 

155:945-959. 

Randall, J. E. 1998. Zoogeography of shore fishes of the 
Indo-Pacific region. Zool. Stud. 37:227-268. 

Richardson, D. L. 1998. Descriptions of the colour patterns of 
juvenile anemonefishes (Pomacentridae, Amphiprion) from 
New South Wales and the Lord Howe - Norfolk Island Region. 
Proc. Linn. Soc. NSW 120:81-86. 

Richardson, D. L. 1999. Correlates of environmental variables 
with patterns in the distribution and abundance of two 
anemonefishes (Pomacentridae, Amphiprion) on an eastern 
Australian sub-tropical reef system. Environ. Biol. Fish 
55:255-263. 

Roberts, C. M., and J. P. Hawkins 1999. Extinction risk in the sea. 

Trends Ecol. Evol. 14:241-246. 
Rousset, F. 2008. Genepop'007: a complete re-implementation of 

the genepop software for Windows and Linux. Mol. Ecol. 

Resour. 8:103-106. 
Santini, S., and G. Polacco. 2006. Finding Nemo, molecular 

phylogeny and evolution of the unusual life style of 

anemonefish. Gene 385:19-27. 
Sea Read. 2009. World's first - hybrid fish at Reef HQ aquarium. 

Issue 26. Available at: www.gbrmpa.gov.au/__ data/.../SeaRead_ 

26_Final_QXPscreen.pdf 
Slatkin, M. 1993. Isolation by distance in equilibrium and 

non-equilibrium populations. Evolution 47:264-279. 
Steadman, D. W. 1995. Prehistoric extinctions of Pacific island 

birds- Biodiversity meets zooarchaeology. Science 

267:1123-1131. 

Storey, J. D. 2002. A direct approach to false discovery rate. J. R. 

Stat. Soc. Ser. B 64:479-498. 
Suchard, M. A., R. E. Weiss, and J. S. Sinsheimer. 2001. Bayesian 

selection of continuous time Markov chain evolutionary 

models. Mol. Biol. Evol. 18:1001-1013. 
Swofford, D. L. 2001. PAUP . Phylogenetic analysis using 

parsimony (and other methods), version 4.0b6. Sinauer 

Associates, Sunderland, MA. 
Tajima, F. 1983. Evolutionary relationship of DNA sequences in 

finite populations. Genetics 105:437-460. 
Tamura, K., J. Dudley, and M. Nei. 2007. MEGA4, molecular 

evolutionary genetics analysis (MEGA) software version 4.0. 

Mol. Biol. Evol. 24:1596-1599. 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1603 



Historic Hybridization between Anemonefish 

Timm, J., M. Figiel, and M. Kochzius 2008. Contrasting patterns 
in species boundaries and evolution of anemonefishes 
(Amphiprioninae, Pomacentridae) in the centre of marine 
biodiversity. Mol. Phylogenet. Evol. 49:268-276. 

van der Meer, M. H., M. Gardner, J. P. Hobbs, and L. van 
Herwerden. 2011. Identification of seventeen microsatellite 
markers for studies on population connectivity in the 
anemonefish Amphiprion mccullochi. Conserv. Genet. Resour. 
DOI: 10.1007/sl2686-011-9517-l 

van Herwerden, L., and P. J. Doherty. 2006. Contrasting genetic 
structures across two hybrid zones of a tropical reef fish, 
Acanthochromis polyacanthus (Bleeker 1855). J. Evol. Biol. 
19:239-252. 

van Herwerden, L., J. H. Choat, C. L. Dudgeon, G. Carlos, S. J. 
Newman, A. Frisch, and M. van Oppen. 2006. Contrasting 
patterns of genetic structure in two species of the coral trout 
Plectropomus (Serranidae) from east and west Australia, 
Introgressive hybridation or ancestral polymorphism. Mol. 
Phylogenet. Evol. 41:420-435. 

van Herwerden, L., J. H. Choat, S. J. Newman, M. Leray, and G. 
Hillersoy. 2009. Complex patterns of population structure and 
recruitment of Plectropomus leopardus (Pisces, Epinephelidae) 
in the Indo- West Pacific, implications for fisheries 
management. Mar. Biol. 156:1595-1607. 



M. H. van der Meer etal. 

Van Oosterhout, C, W. F. Hutchinson, D. P. M. Wills, and P. 
Shipley. 2004. MICRO-CHECKER, software for identifying 
and correcting genotyping errors in microsatellite data. Mol. 
Ecol. Not. 4:535-538. 

Walsh, P. S., D. A. Metzger, and R. Higuchi. 1991. Chelex-100 as a 
medium for simple extraction of DNA for PCR based typing 
from forensic material. Biotechniques 10:506-513. 

Willis, B. L., M. J. H. van Oppen, D. J. Miller, S. V. Vollmer, and 
D. J. Ayre. 2006. The role of hybridization in the evolution of 
reef corals. Annu. Rev. Ecol. Evol. Syst. 37:489-517. 

Winters, K. L., L. van Herwerden, J. H. Choat, and D. R. 
Robertson. 2010. Phylogeography of the Indo-Pacific 
parrotfish Scarus psittacus: isolation generates distinctive 
peripheral populations in two oceans. Mar. Biol. 
157:1679-1691. 

Yaakub, S. M., D. R. Bellwood, L. van Herwerden, and F. M. 
Walsh. 2006. Hybridization in coral reef fishes, introgression 
and bi-directional gene exchange in Thalassoma (family 
Labridae). Mol. Phylogenet. Evol. 40:84-100. 

Zwickl, D. J. 2006. Genetic algorithm approaches for the 
phylogenetic analysis of large biological sequence datasets 
under the maximum likelihood criterion. Ph.D. thesis, The 
University of Texas at Austin, TX. Available at: 
http://repositories.lib.utexas.edu/handle/2152/2666 



1604 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



